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We theoretically study the magnetization inside a normal metal induced in an s-wave 
superconductor/ferromagnetic metal/normal metal/ferromagnetic metal/s-wave superconductor 
{S/El/N/E2/S) Josephson junction. Using the quasiclassical Green’s function method, we show 
that the magnetization becomes finite inside the N. The origin of this magnetization is due to 
odd-frequency spin-triplet Cooper pairs formed by electrons of equal and opposite spins, which 
are induced by the proximity effect in the S/E1/N/E2/S junction. We find that the magnetiza¬ 
tion M{d,0) in the Yean be decomposed into two parts, M{d,0) = M^(d) + M^^(d,0), where 0 
is the superconducting phase difference between the two Ss and d is the thickness of N. The 0 
independent magnetization M^(d) exists generally in S/F junctions, while M^^(d,0) carries all 6 
dependence and represents the fingerprint of the phase coherence between the two Ss in Josephson 
junctions. The 0 dependence thus allows us to control the magnetization in the N by tuning 0 for 
a fixed d. We show that the 6 independent magnetization M^(d) weakly decreases with increasing 
d, while the 0 dependent magnetization d) rapidly decays with d. Moreover, we find that 

the time-averaged magnetization d)) exhibits a discontinuous peak at each resonance DC 

voltage Vn = nhws/^^e (n: integer) when DC voltage V as well as AC voltage Uac(t) with frequency 
ujs are both applied to the S/F1/N/F2/S junction. This is because M^^{d,0) oscillates generally 
in time t (AC magnetization) with dO/dt = 2e[U + Uac(t)]/d and thus (M^^(d, d)) = 0, but can be 
converted into the time-independent DC magnetization for the DC voltage at Vn- We also discuss 
that the magnetization induced in the N can be measurably large in realistic systems. Therefore, 
the measurement of the induced magnetization serves as an alternative way to detect the phase 
coherence between the two Ss in Josephson junctions. Our results also provide a basic concept for 
tunable magnetization in superconducting spintronics devices. 

PACS numbers: 74.45.+c, 72.25.Ba, 74.78.Na 


I. INTRODUCTION 

The proximity effect is an important quantum phe¬ 
nomenon which occurs when a superconductor is at¬ 
tached to non-superconducting materials, where the 
pair amplitude of Cooper pairs in the superconduc¬ 
tor penetrates into the non-superconducting materi¬ 
als [l|. A typical example is the Josephson effect, 
which has been known as one of the macroscopic quan¬ 
tum phenomena, characterized as DC current flow¬ 
ing without a voltage-drop between two superconduc¬ 
tors separated by a thin insulator (7) or normal metal 
m n i. The Josephson critical current in an su¬ 
perconductor/insulator/superconductor or superconduc¬ 
tor/normal metal/superconductor junction forming a 
Josephson junction monotonically decreases with increas¬ 
ing the thickness of / or Y W- 

The proximity effect in 5-wave superconduc¬ 
tor/ferromagnetic metal (S/F) hybrid junctions has 
been extensively studied in the last decade because of 
its fascinating phenomena and potential applications 
to superconducting spintronics (j-Q . Due to the 
proximity effect between S and Y in a S/F junction, 
the spin-singlet Cooper pairs (SSCs) penetrate into 


the F and acquire a finite center-of-mass momentum 
proportional to the exchange splitting between up- and 
down-spin bands in the F. The pair amplitude of SSC 
shows damped oscillation with increasing the thickness 
of F. One interesting phenomena induced by the damped 
oscillatory behavior of the pair amplitude of SSC is a 
TT-state in a S/F/S junction, where the current-phase 
relation in the Josephson junction is shifted by tt from 
that of the ordinary S/I/S or S/N/S junction (called 
0-state) nil. It is expected that the 7r-state can 
be utilized for an element of quantum computing and 
circuit [234^. 

Another intriguing proximity effect in S/F hybrid 
junctions is the emergence of odd-frequency spin-triplet 
Cooper pai rs ( STCs), although the S is an 5-wave super¬ 
conductor [^, [^. Here, the anomalous Green’s func¬ 
tions of spin-triplet components are odd functions with 
respect to the fermion Matsubara frequency oUn- It should 
be noted that the anomalous Green’s functions in bulk 
superconductors are generally even functions with re¬ 
spect to uJn- When the magnetization in the F is uni¬ 
form in a 5/Y junction, not only the SSC, as described 
above, but also the STC composed of opposite spin elec¬ 
trons (i.e., total spin projection on 2: axis being Sz = 0) 
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penetrates into the F due to the proximity effect [H, [29| . 
The penetration length of STC with Sz = 0 (and also 
SSC) into the F is very short and the amplitude of STC 
exhibits a damped oscillatory behavior inside the F with 
increasing the thickne ss of F. T he penetration length is 
determined by = \/J^F>F/h^, which is typically a or¬ 
der of few nanometers Here, and hex are the 

diffusion coefficient and the exchange field in the F, re¬ 
spectively. 

On the contrary, when the magnetization in the F is 
non-uniform in a 5/F junction, the STC formed by elec¬ 
trons of equal spin {\Sz\ = 1) can also be induced in the 
F. This includes cases, for instance, where the F contains 
a magnetic domain wall jsol-l^. the junction consists of 
F multilayers [33-[55|, the interface of S/F junction is 
spin active |56l-l^. and the ferromagnetic resonance oc¬ 
curs Although the pair amplitude of STC with 

\Sz\ = 1 monotonically decreases with increasing the 
thickness of F, the STC with \Sz\ = 1 can propagate into 
the F over a distance of the order of = ^^hDp/^TrkBT 
(T : temperature), which is typically about several dozen 
nanometers [65| . This is approximately 2 orders of mag¬ 
nitude longer than the penetration length of the SSC and 
the STC with = 0. Therefore, the proximity effect of 
STCs with \Sz\ = 1 is called the long-ranged proximity 
effect (LRPE). 

Following the theoretical predictions, the STC in S/F 
hybrid junctions has been confirmed experimentally |6^ 
[ 7 ^ . The obvious way to observe the LRPE induced 
by the STC with \Sz\ = 1 is to directly measure the 
Josephson current in Josephson junctions composed of 
Fs jbbl-lTQj. Indeed, the LRPE has been observed in S/F 
junctions with spin-active interfaces [ggI, [gqI, [70| and in 
S/F multilayer systems with non-collinear magnetization 
alignment between F layers [63, il. Recently, the varia¬ 
tion of superconducting transition temperature (Tc) has 
been observed in F/F1/F2 type spin valve structures as 
the direction of magnetizations in the two ferromagnetic 
metals FI and F2 is changed Ini III. This is also due 
to the LRPE induced by the STC as predicted in the 
previous theoretical calculation j45| . 

An alternative way to prove the STC is to measure 
the spin angular momentum carried by Cooper pairs be¬ 
cause the spin is finite for the STC but is zero for the 
SSC. Several theoretical studies have already addressed 
this issue and examined the magnetization induced by 
the STC in the various geometry of S/F hybrid struc¬ 
tures ( 2 ^, HqI, E, |48| . A F/F/F junction with a spin valve 
structure is a typical geometry of such F/F hybrid struc¬ 
tures. When the magnetizations in the two Fs separated 
by the F are non-collinearlly aligned, not only the STC 
with \Sz\ = 0 but also the STC with \Sz\ = 1 becomes 
finite and induces a finite magnetization inside the F as 
well as the two Fs HilSIil. 

Recently, the magnetization induced by the STC has 
also been studied in Josephson junction type multilayer 
systems, e.g., F/F/F/F, F/F/F/F/F, and rather com¬ 
plex symmetric three terminal F/F/F/F/F/F/F junc¬ 


tions dil. It has been pointed out that such Joseph¬ 
son junctions with metallic ferromagnetic multilayers, es¬ 
pecially, the symmetric three terminal F/F/F/F/F/F/F 
junction may have promising potential for superconduct¬ 
ing spintronics applications with low dissipation ji^ . 
This is because the magnetization in this junction can 
be well controlled by changing the superconducting phase 
difference between the two outmost Fs without Jule heat¬ 
ing. Here, it should be noted that the thickness ds of S 
in the middle layer sandwiched by the two ferromagnetic 
double layers has to be ds Cs /Co (Cs • superconducting 

coherence length) in order to observe clearly the magne¬ 
tization in the middle F layer induced by the STC 
However, in this case, the superconductivity in the mid¬ 
dle F layer is violently suppressed. To prevent this from 
happening, for example, a three terminal Josephson junc¬ 
tion composed of large superconducting electrodes in the 
middle F layer is proposed [jsl, [HU, [73| . 

In this paper, we focus on a much simpler Joseph¬ 
son junction with metallic trilayers, i.e., a F/F1/A/F2/F 
Josephson junction (see Fig. [1]), which nowadays has 
been able to be fabricated experimentally [63, [68 |, [tJ, 
and theoretically examine, by employing the quasiclas- 
sical Green’s function method, the magnetization inside 
the N induced by the odd-frequency STCs composed of 
electrons of equal and opposite spins. Fixing the mag¬ 
netization in F2 along the z direction perpendicular to 
the junction direction {x direction), we show that i) the 
X component of the magnetization in the N is always 
zero, ii) the y component becomes exactly zero when the 
magnetizations in FI and F2 are collinear, and iii) the z 
component is generally finite for any magnetization align¬ 
ment between FI and F2. We also show that the mag¬ 
netization in the N can be decomposed into two parts, 0 
dependent and independent parts, where 0 is the super¬ 
conducting phase difference between the two Fs in the 
F/F1/A^/F2/F junction. The 0 dependent magnetiza¬ 
tion is induced as a result of finite coupling between the 
two Fs, while the 0 independent magnetization always ex¬ 
ists due to the proximity effect in F/F hybrid junctions. 
We find that the 0 independent magnetization decreases 
slowly with increasing the thickness of N^ whereas the 0 
dependent magnetization decays rather rapidly. We also 
investigate the dynamics of the magnetization in the N 
when AC voltage is applied. Because of the AC volt¬ 
age, the superconducting phase difference 0 is now time 
dependent and accordingly the 0 dependent part of the 
magnetization oscillates. However, we find that the 0 
dependent part of the magnetization is converted, when 
it is time averaged, from the oscillating AC character to 
the time-independent DC character at specific DC volt¬ 
ages, depending on the frequency of AC voltage, if DC 
and AC voltages are both applied to the F/F1/A^/F2/F 
junction. Finally, we argue that the magnetization in¬ 
duced inside the N can be large enough to be observed 
experimentally in realistic settings. 

The rest of this paper is organized as follows. In Sec. im 
we introduce a simple F/F1/A^/F2/F junction consisting 
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of metallic trilayers and derive the analytical formulation 
of the magnetization induced inside the N on the basis of 
Usadel equation in the diffusive transport limit. It is clear 
from this analytical formulation that the magnetization 
in the N is indeed induced by the odd-frequency STCs. 
In Sec. imi we show the results of the magnetization as 
a function of the thickness of N for different magneti¬ 
zation alignments of the two Fs. The 0 dependence of 
magnetization, including the dynamics when AC voltage 
is applied, is also discussed. Finally, the magnetization 
induced by the STCs is estimated for a typical set of re¬ 
alistic parameters in Sec. [IVl The summary of this paper 
is given in Sec. O The spatial dependence of anomalous 
Green’s functions in the N is discussed in Appendix [Al 
and the local magnetization density induced inside the 
N is examined in Appendix [B1 

II. JUNCTION AND FORMULATION 

After introducing the Josephson junction studied, we 
first formulate for this junction the anomalous Green’s 
functions in the diffusive transport limit on the basis of 
the quasiclassical Green’s function method and then de¬ 
rive the analytical formulae of the magnetization induced 
inside the N. 


A. S/F1/N/F2/S junction 

As depicted in Fig. (TJ we consider the S / FI/ N/ F2/S 
junction made of normal metal {N) sandwiched by two 
layers of ferromagnetic metal (FI and F2) attached to s- 
wave superconductors (Fs). We assume that the magne¬ 
tization in F2 is fixed along the z direction perpendicular 
to the junction direction {x direction), while the FI is a 
free layer in which the magnetization can be controlled 
by an external magnetic field, pointing any direction in 
the yz plane, parallel to the interfaces, with (p being the 
polar angle of the magnetization. We also assume that 
the magnetizations in FI and F2 are both uniform. The 
thicknesses of F, FI, F2, and N are ds, dpi, dF 2 , and 
d, respectively, with L = d F dpi, Fp = F dF 2 , and 
Fs = Fp -1- ds- Furthermore, we assume that ds is much 
larger than the superconducting coherent length ^s- 


B. Anomalous Green’s functions 

In the diffusive transport region, the magnetization in¬ 
side the N is evaluated by solving the linearized Usadel 
equation in each region m (= FI, A, and F2) 


ihDmdlFix) -i‘2h\L0n\ Fix) +2A{x)-sgn{LOn)hl^{x) I^Ty,F{x)'^-sgn{ujn)K^{x) %,F{x) = 6 , ( 1 ) 


where Dm is the diffusion coefficient in region ujn = 
(2n + l)T:k^T/h with n = 0,±1,±2 ,'-' is the fermion 
Matsubara frequency, sgn(A) = A/|A|, and fy(^z) is the y 
(z) component of Pauli matrix. We assume that diffusion 
coefficients in FI and F2 are the same, i.e. Fpi = Ff 2 = 
Fp. Note also that {Q, R} = QR + RQ, [Q, R] = QR — 
RQ^ and 0 is null matrix. The anomalous part f'^ of the 
(2 X 2) quasiclassical Green’s function [56| is given by 


i.e., 


A(a;) = 



, -ds < X < 0 

, Fp < X < Fs 

other 


(3) 


r {x) 


( /-ft (^) fn i 

( -ftS (x) + ifl^ {x) {x) + /j™ {x) Y . 

'v -/r (x)+m (x) m (^)+ if^y i^) 


where the ujn dependence is implicitly assumed. Notice 
that f^{x) is the anomalous Green’s function for the 
SSG, whereas and fl^{x) represent the anoma¬ 

lous Green’s functions for the STC with \Sz\ = 1 and 
\Sz\ = 0, respectively. The s-wave superconducting gap 
A{x) is finite only in the F and assume to be constant. 


The exchange field he^{x) = {h^^{x), hl^{x), h^^{x)) due 
to the ferromagnetic magnetization in the Fs is described 
by 


dex (^) 


KF + KJz, 

0 , 


0 < X < dpi 

L < X < Fp, 
other 


(4) 


where = dexi sin(/9, = dexi cosip (see Fig. [T]), and 

ey(^z) is a unit vector in the y (z) direction. We assume 
that dexi and dex 2 are both positive. 
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FIG. 1: (Color online) Schematic illustration of the 

S/Fl/N/F2/S junction studied, where the normal metal (N) 
is sandwiched by two ferromagnetic metals (FI and F2) at¬ 
tached to s-wave superconductors (Ss). Arrows in FI and F2 
indicate the direction of ferromagnetic magnetizations. While 
the magnetization in F2 is fixed along the z direction, the FI 
is assumed to be a free layer in which the magnetization can 
be controlled by an external magnetic field within the yz plane 
with if being the polar angle of the magnetization, ds, dpi, 
dF 2 , and d are the thicknesses of S', FI, F2, and N, respec¬ 
tively, with L = d F dpi, Lp = L + dp 2 , and Ls = Fp + ds. 
We assume that the magnetizations are uniform in both FI 
and F2 layers, and that ds ^ ^s- 


To obtain the solutions of Eq. o , we impose appro- 
priate boundary conditions [75|, i.e., 



) = Fix) 

x=0 

x=0 

Fix) 

: 

y = F(^) 

C = «F1 

x^d-pi 

Fix] 

1 ^ = Fix) 

x=L 

5 

x=L 

Fix) 

, = Fix) 

X = Lf ■ 

5 

x=Lf 

dxFix) 

= —dxf 

C = dFl 7P 

'^ix) 


( 5 ) 

(6) 

( 7 ) 

(8) 
( 9 ) 


J 


and 


Fda;F{x) 

7f 


x=L 


= d,f\x) 


x=L 


( 10 ) 


where yp = cr-p/a-^ and crF(N) is the conductivity of FI 
and F2 {N). Moreover, in the present calculation, we 
adopt the rigid boundary condition 


<TS CS ’ 


( 11 ) 


where as is the conductivity of S in the normal state and 
Cfi( 2 ) = i//iT)F//iexi( 2 ) 0- Assumiug that ds > Cs, the 
anomalous Green’s function in the Ss attached to FI and 
F2 can be approximately given as 


/f(a;)U=0(LF) = -Ty 


^L(R) 


+ I^L(R)P ' 


( 12 ) 


where Al(r) = Ae^^^cR) (/\. ^eal) and dL(R) is the su¬ 
perconducting phase in the left (right) side of Ss (see 
Fig. [T]). 

Assuming that d^i/^Fi ^ 1, we can preform the Tay- 
lar expansion for f^^{x) as follows IHIzi: 


f^{x) « fHdFi)+ix-dFi)d,f^x) 


X=dFi 


+ 


{x dpi) o2 Ffi 




X = dFl 


(13) 


Using the boundary conditions given in Eqs. (|5]) and (|9]) 
for Eq. (p!3|) and substituting Eq. ([13]) into Eq. ([T]), f^^{x) 
can be approximately expressed as 


f\x) 


— dj^ix) 

7f 




- zsgn(w„) 


(x - dpif 
2dFp 


(*L{'7yd®^(0)} + f^,/®^(0) ) 


(14) 


Here we also assume that the exchange field dgxi in ary conditions given in Eqs. (|H|) and (flO]) . f^‘^{x) can be 
the FI is much larger than /crT and thus the term approximately expressed as 
d\ujn\f^^{x) is neglected in Eq. ([T4|) . 

Similarly, assuming that dp 2 /Cf 2 ^ 1, we can perform 
the Taylor expansion for f^‘^{x) and, using the bound- 
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FHx) « dj^ix) + PHlf) + ^ 

7f x=l 7f 


. ^ex2^ 


F 2 


d.f M__^ + »g„K) 


,/“(ip) 


- zsgn(w„) 


{x Z/) /lex2 
2fiDF 


Tz,P'^{Lf) 


(15) 


where hex 2 ^ feT is also assumed. 

The general solutions of f^{x) are given as 

f^{x) = 

(16) 

and /cn = ^/2\uJn\/D^• Applying the boundary condi¬ 
tions given in Eqs. ([6]) and ([fj) to Eqs. (pT|) and (pFj). and 
also using the results in Eqs. (p!T|) and (pF]). we can obtain 
the anomalous Green’s functions in the N as 

and 




(17) 


where 



f^{x) = /f(x)±/,7x) 

(18) 



J 


/. (^) = 


+ i 


' sinh[fcN(a; — L)] — cosh[fcN(a; — L)] K^^{d) 

7f 

sinh[fcN(a; - dpi)] + cosh[k^{x - dpi)] Ku^{d), 




■ Ar 




f^y{x) = Sgn(w„) 


Al fe|x<^Fl 
HDy 


7f 


sinh[/cN(^ — L)] — 


feN<^F2 

7f 


cosh[/cN(a:^ — L)] 


K.Ad). 


(19) 

( 20 ) 


and 


ffAx) 


= sgn(w„) 
- sgn(w„) 


Al /tex4l 

Eu„ dDp 
Ar /lex2C^F2 
Eu,„ hD-p 


sinh[fcN(a; — L)] — cosh[^ n( a: — L)] K^^{d) 

7f J 

sinh[fcN(a: - ^Fi)] + cosh.[kpi{x - dpi)] Ku„{d), 

7f j 


( 21 ) 


where 

E^^ = v/(Sa;„)2 + A2 (22) 

and 

d^uj^W = -+- cosh(A)Nd) 

V 7f 7f / 

+ fl + ^^)sinh(fcNd). (23) 

V If If J 

Prom Eqs. (P1?| - (PT|) . it is immediately found that that 
ff{x) describing the SSC is an even function with re¬ 
spect to ujni whereas /^^ 7 ^)(^) describing the STC is an 
odd function with respect to ujn since is pro¬ 

portional to sgn(co’^). Hence, f^y{tz)^^) represents the 
odd-frequency STC. 


It should be emphasized here that 

lim /t7^)=0 (24) 

/lexl^O ^ 

and 

, lim ffz{x)=Q, (25) 

^exl ,^ex2^U 

whereas ff{x) is generally finite independently of hexi 
and hex 2 - This is due to the fact that the presence of 
F layers are essential to induce the STC [^. On the 
contrary, the SSC is always induced inside the N in S/N 
junctions as well as more complex S'/A^/F junctions 
Notice also that i) f^{x) = 0 because the exchange field 
in the FI does not have the x component, and ii) f^{x) 
is exactly zero when (^ = 0 or tt as f^{x) cc The 
spacial dependence of anomalous Green’s functions in the 
N is discussed in Appendix [Al 
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C. Induced magnetization in normal metal 

Within the quasiclassical Green’s function method, the 
magnetization M(d, 0) induced inside the N is given 
as 


M{d,e) = {M^{d,0),My{d,e),M,{d,e)) 

A 

= — m{x,0)dx, (26) 

^ JdFi 

where 0 = — is the superconducting phase difference 

between the outmost ^Ss in the junction and 


i{x,0) = {mx{x,0),my{x,0),mz{x,0)) 
= -giJ.BT^NFkBT'^sgn{LOn)lm 


(27) 


with 


= (/te(a;),-/ty(a;),/('^(a;)) 


(28) 


Here, m(x, 0) is the local magnetization density in the N, 
g is the g factor of electron, fiB is the Bohr magneton, 
and A and V = Ad are the cross-section area of junction 
and the volume of N, respectively. In the quasiclassical 


J 


Green’s function method, the density of states N-p per 
unit volume and per electron spin at the Fermi energy is 
assumed to be approximately the same for up and down 
electrons in the N liS-Iiil. 

It is apparent in Eq. (|27ll that f^{x) and ft^(x) are 
both required to be nonzero to induce finite rh{x,0). 
However, as described in Sec. imi nonzero ft^{x) oc¬ 
curs only when F layers are involved in the junction and 
ft^{x) = 0 whenever ff{x) = 0 for Al = Ar = 0 . 
Therefore, the origin of the magnetization in the N is 
considered to be due to the STGs induced by the proxim¬ 
ity effect [ill, liO, Eil . Note also that because of = 0 

(see Sec. imi), mx{x,0) and thus Mx{d,0) are always 
zero. Therefore, in the following, we only consider the 
y and z components of M{d^0). More details of m{x,0) 
are examined in Appendix [0 

Substituting Eqs into Eq. (|27|) and perform¬ 

ing the integration with respect to x in Eq. (|26|) . we can 
obtain the y and z components of the magnetization in¬ 
duced inside the N. The y component My{d,0) of the 
magnetization is decomposed into two parts. 


My{d,0) = Mlid)+Ml\d,d), (29) 


where 


Mhd) = -guBkBT- 


ttTVfA^ hyjl,^KUd)Fl(d) 


2d HDb 


E- 


UJri V / ' 


(30) 


and 


Ml\d,0) = 5MB 


ttTVfA^ hyjl, ^ Ky(d)Fy (d) 


2d HDb 






COS^. 


(31) 


Here, we have introduced 


[1 _ cos(2fcNd)] 
7f 


^ / ^N^F2 Y 

V 7F y 


- 

2 


/ h^d^\ ^ 
V 7f y 


sinh(2/cNd) 


(32) 


and 


(^) = ( ^ + feN^Fl feN^F2 j 
7f 7f 


^ _ ^Ndpi k]sidB2 _ f ^N^Fl ^N^F2 \ 
7f 7f 

I- 


H—j 
7f 7f y . 


smh{k^d). (33) 


Similarly, the z component Mz (d, 0) of the magnetization where 
is decomposed into two parts. 


Mz{d,0) = Ml{d) F Ml\d,0), 


(34) 


Ml{d) = ggBkBT 


ttNbA^ ^ K^Jd) 
2d ^ 


^ex^Fl oa / I ^ex2d|2 


HDt^ 


■KAd)F 


HDt^ 


KAd) 


(35) 
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and 




( 36 ) 


Here, we have also introduced 

R£>(<I) = + 1 - 


^NC^F1(2) Y 

7f J 


k]^d-pT(o') , ^ 1 

A:Nd - ^^ cosh(2feNd) - - 

7f 2 


kfjd] 


’fi(2) Y 


7f / 


sinh(2/cNd). (37) 


r 


The 6> independent part of the magnetization, i.e., 
Mlid) and Ml{d), is due to the proximity effect common 
in 5'/F junctions, similar to the one inducin g t he STCs 
in F/S/F and S/F/F junctions [HI, |40|, EH, EH, Ell, E3 • 
On the other hand, the 0 dependent part My^{d,0) and 
Ml^{d,0) of the magnetization is induced by the cou¬ 
pling between the two Ss in the junction. Therefore, 
^y\z) only when ferromagnetic metallic 

multilayers constitute the Josephson junction [48|. It 
should also be noticed that My{d,0) becomes zero when 
= 0 or TT since My{d) and Mii(d, 6>) are both propor¬ 
tional to = hexi sin(/:?. In contrast, Mz{d^ 0) is gener¬ 
ally nonzero for any (p. 


III. RESULTS 

A. Thickness dependence of magnetization in 
normal metal 

Let us first numerically evaluate My{d, 0) and Mz{d^ 0) 
in the N obtained in Eqs. (|29|) - (|37|) . For this purpose, 
the temperature dependence of A is assumed as 

A = Aotanh ^1.74^^ - 1^ , (38) 

where Aq is the superconducting gap at zero tempera¬ 
ture and Tq is the superconducting transition temper¬ 
ature 0- Figures [2]-0] show the typical results of the 
magnetization in the A as a function of thickness d of 
the N normalized by = ^/hD^^J2^^k^Fc. 

Figure [2] represents the results for (^ = 0 where the 
magnetizations between FI and F2 are parallel. As 
shown in Fig. [2{a), the y component My{d^0) of the 
magnetization is exactly zero since f^{x) contributing 
to My{d,0) is zero in the parallel magnetization config¬ 
uration. On the other hand, the 2 ) component Mz{d^0) 
has a finite value, as shown in Fig. Efb), because ffz{^) 
contributing to Mz (d, 0) is non-zero in the parallel mag¬ 
netization configuration. Furthermore, the induced mag¬ 
netization Mz{d^ 0) is found to be negative, i.e., opposite 
to the magnetizations in FI and F2. It is also found 
in Fig.[2jb) that \Mz{d^0)\ monotonically decreases with 
increasing d for d > but M](d) decays rather slowly 



FIG. 2: (Color online) (a) The y component My(d,0) and 
(b) the z component Mz{d,0) of the magnetization in the 
N for 9 ? = 0, corresponding to the parallel magnetization 
configuration between FI and F2. For other parameters, we 
set TjTc = 0.3, ^ = 0, 7f = 0.1, dpi/^n = 0.3, dF 2 /^D = 0.2, 
dexi/Ao = 30, and /iex 2 /Ao = 20. For comparison, Ml{d) 
and Ml^(d,0) are also plotted separately in (b). 


as compared with M]^(d, d). The difference of the decay 
rates for d ^ as well as the small d behavior in M](d) 
and Mj^(d, d) will be further discussed below. 

Figure [3] shows the results for p = ir where the magne¬ 
tizations between FI and F2 are antiparallel. As shown 
in Fig. [3l^a), the y component My{d^0) of the magne¬ 
tization is still exactly zero since f^{x) contributing to 
My{d,0) is zero also in the antiparallel magnetization 
configuration. However, the 2 ) component Mz{d^0) is fi¬ 
nite and decreases monotonically with increasing d for 
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FIG. 3: (Color online) (a) The y component My{d,0) and 
(b) the z component Mz{d,0) of the magnetization in the N 
for (p = TT, corresponding to the antiparallel magnetization 
conhguration between FI and F2. Other parameters are the 
same as in Fig. [2] For comparison, Ml{d) and Ml^{d,0) are 
also plotted separately in (b). 



FIG. 4: (Color online) (a) The y component My{d,0) and 
(b) the z component Mz{d,0) of the magnetization in the N 
for p = 7r/2, corresponding to the case where the magnetiza¬ 
tion in FI is perpendicular to that in F2. Other parameters 
are the same as in Fig. O For comparison, M^(^)(d) and 
^y\z) (<^5 plotted separately. 


d > [see Fig. [3l[b)]. It is also noticed in Fig.[3][b) that 
the sign of Mz{d^ 0) is positive and is opposite to the one 
for (/? = 0 [Fig.[2][b)]. The sign reversal of Mz{d, 0) will be 
further discussed below. It is also observed in Fig. [3][b) 
that Ml{d) decays slowly with respect to d as compared 
with M^^((i, 6>), similarly to the case when the magnetiza¬ 
tions in FI and F2 are parallel. It should be noted how¬ 
ever that Mz {d^ 0) becomes exactly zero in the antiparal¬ 
lel magnetization configuration when dpi = dY 2 and \hl^\ 
= |^ex 2 | because in this case M^(d) and M^J^{d^0) are 
both zero, as seen in Eqs. dSSli-dSZl) [see also Eqs. m- 
(031)]. 

EigurelH shows the results for p = n/2 where the mag¬ 
netization in FI is perpendicular to that in F2. As shown 
in Eig.|ll[a), the y component My{d, 0) of the magnetiza¬ 
tion is now finite since ffy{x) contributing to My{d, 0) is 
nonzero in this case. It is also found in Eig. |4]that both 
\My{d^ 0) \ and \Mz{d, 0)\ decrease monotonically with in¬ 
creasing d for d > ^D- Moreover, it is clearly observed 
that the decay rate of |M^^^^(d)| with respect to d is 

slower than that of ^)|. This is similar to the 

other cases discussed above in Eigs. [2]and[3l 

It is now instructive to consider limiting cases for the 


magnetization M^^^^(d) and 6>) induced inside 

the N and analyze the qualitative behavior with respect 
to the thickness d of the N. Eor T ^ Tq and d ^ the 
y components My{d) and M^{d^6) of the magnetization 
are approximately given as 



Mlid) 


(39) 


27rkBT HDb d 

and 

Ml\d,e) ^ 

tt/cbT HDb 

(40) 

whereas the 2 ; components Ml{d) and Ml^{d^0) 
magnetization are approximately 

of the 


Mlid) 

_ ^/iBApA^ Fex 

271 k^T hF)^ d 

(41) 

and 

Ml\d,e) « 

f“ CO. 

TT k^ T HDy 

(42) 

where 

He. 

^ex^Fl Aex2^P2 

(43) 
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and = \/hDj<\ /27r T at T ~ Tc). It is 

immediately found in Eqs. (|39])-(j42j) that de¬ 

creases rather slowly, i.e., algebraically, as 1/d, whereas 
^y\z) decays exponentially. This is indeed compa¬ 
rable with the numerical results shown in Figs. [2]-|4l 

Next, we discuss in the same limiting case the sign 
change of Mz (d, 0) by flipping the magnetization direc¬ 
tion from the parallel to the antiparallel configuration in 
FI and F2. For this purpose, we focus on idex appear¬ 
ing in Eqs. (|4T]) and (|42j), the definition being given in 
Eq. (|43|) . In the case of parallel magnetization configu¬ 
ration, Mz{d,0) is always negative, as shown in Fig. [2l 
simply because is positive (assuming that \0\ < 7r/2). 
On the other hand, in the case of antiparallel magnetiza¬ 
tion configuration, is negative. Therefore, Mz{d^0) 
becomes positive when is larger than dex 2 d |2 5 

shown in Fig. [3l 

Let us now consider the opposite limit, i.e., d <C 
at T ^ Tc. In this limit, the y components M^(d) and 
My^{d,0) of the magnetization are approximately given 
as 


Mlid) ^ 


X 


FI f dpi 

tt/cbT HDf V^n7f 
dF2 f dF2 ^ d \ 

Cf7f VCn7f Cn/ 


dF2 \ 

Cn7f/ 

(44) 


J 


Ml\d,9) 


giiBNpA"^ / ^ex2dp2 \ / dpi 

tt/cbT \ HDf dFp / VCn7f 


It is therefore readily noticed in Eqs. (|44|) - (|Tf|) that 
M^(^)(d) and M^|^^(d, 0) are linearly dependent on d and 
their slopes are determined by the signs of and 
This is in good qualitative agreement with the numerical 
results shown in Figs. [2]-[4l 


B. 6 dependence of magnetization in normal metal 

In the previous section, we have focused on the d de¬ 
pendence of the magnetization induced inside the N. 
Here, we shall demonstrate that the magnetization can 
also be controlled by the superconducting phase differ¬ 
ence 0 in the two Ss. The most simplest way to tune 0 ex¬ 
perimentally is to apply DC bias current to the junction, 
in which the DC Josephson effect can be detected [t^ . 

Figure [5] shows the y and 2 ; components My{d^0) and 
Mz (d, 0) of the magnetization induced inside the as a 
function of 0 for different values of yp. Figure [5] clearly 
demonstrates that the magnetization can indeed be con¬ 
trolled by tuning 0. It should also be noticed that the 


and 


Ml\d,0) 


X 


^L^pi f <^Fi ^ dF2 \ 
tt/cbT HDf V^n7f Cn7f / 

dpi dF2 / dpi ^ dF2 \ d 
_Cn7f Cn7f V2^n7f 2^n7f/ ^n_ 

(45) 


whereas the 2 ; components Mj(d) and M]^(d, d) of the 
magnetization are approximately 


Mlid) 


and 


^/iB^F^^ dex2C^; 


FI 


di 




nksT hDp 


f d^i 
VCnTf 


k 


(. 


'F2 


+ — 
^ex2 V^N7F Cn 


<^F2 \ 

Cn7f/ 


-2 


dpi d] 


(. 


d 


/ dF2 Y 
7n7f/ 




(46) 


di 


, ■ )- 

Cn7f / LCn7f Cn7f V^^nTf 2^n7f/Cn, 


id. 


(47) 


magnitude of the magnetization increases with decreas¬ 
ing 7 f. Therefore, yp is an important parameter to in¬ 
crease the magnetization induced inside the N. 


C. Dynamics of magnetization in normal metal 

Next, let us discuss an alternative way to control the 
magnetization induced inside the N. Here, we consider 
the 5'/Fl/A^/F2/5'junction subject to both DC and AC 
external fields, described by the voltage bias model [78|, 
as schematically shown in Fig.[6l In this case, the super¬ 
conducting phase difference d evolves in time t according 
to the following well known formula: 

9(t) = 9o + sm{ujst), (48) 

tl tlUJs 

where do is a time independent constant, V is the DC 
voltage, and vs and ujs are the amplitude and frequency 
of the AC voltage, respectively. 

Substituting Eq. (|48|) into Eqs. (|3T]) and (|3^ . and using 
the generating function of Bessel functions, we can easily 






































LO 



FIG. 5: (Color online) The 0 dependence of the magnetizations, (a) My{d,0) and (b) Mz{d,0), induced inside the N for 
(f = 7r/2, corresponding to the case where the magnetization in FI is perpendicular to that in F2. We set d/^r> = 1 for three 
different values of 7 f indicated in the figures. Other parameters are the same as in Fig. [2l 


find that the 0 dependent parts of the magnetization in 
the N are given as 


FI F2 



^h“© 

V v^cos{oj^t) 


FIG. 6: (Color online) Schematic set up of the 

S/Fl/N/F2/S junction to observe dynamics of the magne¬ 
tization induced inside the N. This set up is based on the 
voltage bias model [t^, where V is the DC voltage, and vs 
and Co’s are the amplitude and frequency of the AC voltage, 
respectively. 


J 


and 


Ml\d,V,t) 


gjj^BkBT 


hy^dl^ 
2d hF)^ 


Ml\d,V,t) 


—gidEksT 


2d 


I k 

V hOv 


knEl^ 




(49) 


(50) 


sides. In the above, we have also introduced 


T{V,t) = 


E 


/ 2ttv^ \ 

V / 


for the y and 2 ; components, respectively, where the V 
and t dependence is explicitly shown in the left hand 


m= —00 


cos [00 + (wj - mu}s)t ], 

( 51 ) 
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where uj = 2eV/h is the Josephson frequency, = h,/2e 
is the flux quantum, and Jm{x) is the Bessel function of 
the first kind (m: integer). 

Let us now consider the time averaged quantity 

SMll^{d,V) = \im^^ r dt (52) 

^ 0 

As shown in Eqs. (|^ and (|5Q]) . the 0 dependent part of 
the magnetization clearly oscillates in t. Therefore, the 
time averaged magnetization V) is zero except 

for ujj = mujs- This is simply because the coefficient pro¬ 
portional to t in the cosine function of Eq. m becomes 
zero only when this condition is satisfied. Indeed, the 
characteristic feature of E) is clearly found in 

Fig.El i.e., V) showing non-zero values only at 

V = mfico’s/2e. 

IV. DISCUSSION 

Einally, we shall approximately estimate the amplitude 
of the magnetization induced inside the N. As shown in 
Eigs. [SHU the magnetization in the N has a finite value 
in the length scale of ^d- In dirty normal metals, is in 
a range of several dozen to several hundred nanometers 
since is about 0.01-0.1 m^/s and Tq is assumed to be 
a few kelvins [65|. As indicated in Eigs. [2]-[5l andO the 
amplitude of the magnetization is estimated to be one 
to two orders smaller than Mq = ^/iBA^FfeTc^D/'?F 2 ' 
When we use a typical set of parameters, i.e., the density 
of states at the Eermi energy 

^ ^-4 X 10^^ eV“^m“^ (53) 

with the Eermi energy e y ^ 5 eV [s^ (m: the electron 
massb = 9 K for Nb [65|, = 100 nm, and ^f 2 = 5 

nm p, we can estimate that Mq is approximately 
31000 A/m. It is therefore expected that the magneti¬ 
zation induced inside the N can be detected with the 
magnetization measurement by SQUID j^ . 

V. SUMMARY 

We have calculated the magnetization inside the N 
in the S /FI/N / F2/ S Josephson junction based on the 
the quasiclassical Green’s function method in the diffu¬ 
sive transport limit. By solving the Usadel equation, we 
have found that finite magnetization is induced inside the 
N. We have shown that the magnetization is due to the 
odd-frequency STCs formed by electrons of equal and op¬ 
posite spins, which are induced by the proximity effect 
in the S/ FI/N/F2 /S junction. Eixing the magnetization 
in F2 along the 2 ) direction perpendicular to the junc¬ 
tion direction {x direction), we have shown that i) the x 
component of the magnetization in the N is always zero. 


ii) the y component is exactly zero when the magneti¬ 
zation direction between FI and F2 is collinear, and iii) 
the 2 ) component is generally finite for any magnetization 
direction between FI and F2. 

Decomposing the induced magnetization into 0 inde¬ 
pendent and dependent parts, we have found that the 
0 independent part of the magnetization decays slowly 
with increasing the thickness of the V, whereas the 0 de¬ 
pendent part of the magnetization decays rather rapidly. 
While the 0 independent part of the magnetization is 
generally induced even in the S/F junctions due to the 
proximity effect, the 0 dependent part of the magneti¬ 
zation results from the finite coupling between the two 
Fs in the S/ FI/ N/ F2/ S Josephson junction. We have 
also found that the time averaged magnetization in the N 
exhibits discontinuous peaks at particular values of DC 
voltage when DC and AC voltages are both applied to the 
S/FI/N/F2/S junction, implying that the AC magne¬ 
tization oscillation can be converted into the DC compo¬ 
nent. We have discussed that the magnetization induced 
inside the N can be large enough to be observed in typ¬ 
ical experimental settings. It is therefore expected that 
a Josephson junction composed of ferromagnetic metal¬ 
lic multilayers such as the one studied here can have a 
promising potential for low Joule heating spintronics de¬ 
vices, where the magnetization can be controlled by vary¬ 
ing the superconducting phase difference d. 
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Appendix A: Spatial dependence of anomalous 
Green’s functions inside normal metal 

In this Appendix, we shall discuss the spatial depen¬ 
dence of the anomalous Green’s functions inside the N. 
The analytical solutions are obtained by solving the lin¬ 
earized Usadel equation (see Sec. Ill Bp and are given in 
Eqs. ([T^ - (l2T]) . Eigure [8] shows the typical results of 
the anomalous Green’s junctions inside the N for three 
different magnetization alignment between FI and F2, 
parametrized by (see Eig. [1]). As shown in Eig. [8fa), 
f^{x) does not depend on cp and exhibits symmetric be¬ 
havior with respect to x about the center of the N. The 
(f independence is simply because ff{x) represents the 
SSC which can be induced even without F layers [see 
Eq. 1^]. 

Eigure [8][b) shows the spatial dependence of the 
anomalous Green’s function f/^{x), corresponding to the 
STC with \Sz \ = 1. Eor the collinear magnetization align¬ 
ment, i.e., (/:? = 0 or TT, f/^{x) is exactly zero because 
f/^{x) is proportional to the y component of the magne- 
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FIG. 7: (Color online) Time average, (a) SMy{d,V) and (b) SMz{d,V), for the 0 dependent part of the magnetization 
Ml^{d,V,t) and Ml^{d,V,t) induced inside the N ior (p = 7r/2, corresponding to the case where the magnetization in FI is 
perpendicular to that in F2. We set d/^u = 1, = 0, and 7ri’s/^o^<t;s = 1- Other parameters are the same as in Fig. [21 


tization in FI [see Eq. (|2Q|) ]. Therefore, the local mag¬ 
netization density my{x^0) and thus the magnetization 
My{d^ 0) in the N is exactly zero in this case (see Figs. [2] 
and|3|). In contrast, /^^(x), corresponding to the STC 
with | 5 '; 2 | = 0, is generally finite, as shown in Fig. IHfc). 
Therefore, the magnetization Mz{d,0) inside the N is 
generally finite, as shown in Figs. |2]-|4l 

It should be noted here that although the analytical 
solutions in Eqs. indicate their exponential de¬ 

pendence with respect to x, Eigs. IHl[b) and |Hl[c) suggest 
that f^{x) for (/? 7 ^ 0 , 7 r and almost linearly. 

This seemingly linear dependence is simply because of the 
parameter set chosen in Eig.jSl where /cnCd = ^T/Tc ~ 
0.5 and hence k-^{L — x) and k-^ix — d-pi) in the exponents 
are no larger than 0 . 6 . 

Appendix B: Local magnetization density inside 
normal metal 

In this Appendix, we will first provide the analytical 
form of the local magnetization density induced inside the 
N and examine the ip dependence. Within the quasiclas- 
sical Green’s function method, the local magnetization 
density m{x^0) inside the N is obtained by substituting 
Eqs. (fT^ - ([ 2 T]) into Eq. (EZl). The y component my{x, 6^) 
of the local magnetization density can be decomposed 
into 0 independent and dependent parts 

my{x,0) = m}y{x) + rpj{x,0), (Bl) 


where 

j^y y2 _ a 2 

m\{x) = F. 

iOn 

(B2) 

and 

hv d‘^ 

m^px,d) = -ggB'xkBTNF-^^^ 

X ix)Ru,^ (x) cos 0{E3) 

Here, we have introduced 

F^^{x) = smh[k-^{x —— cosh[feN( 3 :^ —E)], (B4) 

7f 

and 

R^^{x) = sinh[/cN(x — d-pi)] + cosh[/cN(x — d^i)]- 

7f 

(B5) 

Similarly, the 2 ; component mz{x,0) of the local magne¬ 
tization density can be decomposed into two parts 

mz{x, 0) = m^(x) -h 0), (B 6 ) 


where 
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FIG. 8: (Color online) Spatial dependence of the anomalous 
Green’s functions, (a) fflx), (b) f^y{x), and (c) inside 

the iV at cjn = TvkBT/h for (p = 0, 7r/2, and tt indicated in (a). 
We set the thickness d of the N to be and = 0. 

Other parameters are the same as in Fig. [2] Note that the N 
layer is located in 0.3 < x/^b < 1.3 for this parameter set. 


ml{x) = -gjiBT^kBTN-F 

flJjY 


A2 




{x) - gfiET^ksTN f 


J2—KlJd)Rl^ix) 


(B7) 


and 
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FIG. 9: (Color online) (a) The y component my{x, 0) and (b) 
the z component mz{x,0) of the local magnetization density 
in the N for 9 ? = 0, corresponding to the parallel magnetiza¬ 
tion configuration between FI and F2. We set the thickness 
d of the N to be and other parameters are the same as in 
Fig. [2 For comparison, m\{x) and ni}z{x,0) are also plotted 
separately in (b). 


Figure [ 9 ] shows the numerical results for (^ = 0 where 
the magnetizations between FI and F2 are parallel. As 
shown in Fig. [9](a), the y component my{x^0) of the lo¬ 
cal magnetization density is exactly zero because ffy{x) 
contributing to my{x, 0) is zero in the parallel magnetiza¬ 
tion configuration. On the other hand, the 2 ; component 
mz{x^0) of the local magnetization density has a finite 
value, as shown in Fig. (H^b), since contributing 

to mz{x, 0) is nonzero in the parallel magnetization con¬ 
figuration. Furthermore, the induced local magnetization 
density mz{x^ 0) is found to be negative, i.e., pointing the 
opposite direction to the magnetizations in FI and F2. 
It should also be noticed that both m\{x) and m}^{x^0) 
exhibit generally nonmonotonic behavior with respect to 

X. 

Figure [To] shows the numerical results for (p = ir where 
the magnetizations between FI and F2 are antiparallel. 



X / 4d 


FIG. 10: (Golor online) (a) The y component my{x,0) and 
(b) the z component mz{x^0) of the local magnetization den¬ 
sity in the A for (/? = tt, corresponding to the antiparallel 
magnetization configuration between FI and F2. Other pa¬ 
rameters are the same as in Fig. (9} For comparison, m\{x) 
and ctre also plotted separately in (b). 


As shown in Fig. cni: a), the y component my{x^ 0) of the 
local magnetization density is exactly zero since ffy{x) 
contributing to my{x^0) is zero also in the antiparal¬ 
lel magnetization configuration. On the other hand, as 
shown in Fig. fTOTb). the 2 : component mz{x, 0) of the lo¬ 
cal magnetization density has a finite value since //^(x) 
contributing to mz{x,0) is nonzero in the antiparallel 
magnetization configuration. Furthermore, as opposed 
to the case for (/? = 0, the induced local magnetization 
density mz{x,0) changes the sign from positive to neg¬ 
ative with increasing x. Note also that mY{x^0) is ex¬ 
actly zero for the special case when dpi =^iF 2 and \h^^\ 
= |^ex 2 |, as shown in Fig. [TTl and thus the local magne¬ 
tization density is no longer dependent on 0. 

Finally, Fig. [12] shows the results for (p = 7r/2 where 
the magnetization in FI is perpendicular to that in F2. 
As shown in Fig. da), the y component my{x,0) of 
the local magnetization density is now finite because 
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f^{x) contributing to my{x^0) is nonzero in this case 
[see Fig. [Hl[b)]. Similarly to the previous cases for ip = 0 
and TT, the 2 ; component mz{x^ 0) of the local magnetiza¬ 
tion density is also finite [Fig. fT^ b)]. 


FIG. 11: (Color online) The z component mz{x,0) of the 
local magnetization density in the N for 9 ? = tt, corre¬ 
sponding to the antiparallel magnetization configuration be¬ 
tween FI and F2. We set dpi/fo = dF 2 /fo = 0.2 and 
/lexi/Ao = /iex 2 /Ao = 20. Other parameters are the same 
as in Fig. [9l For comparison, m\(x) and ni}z(x,0) are also 
plotted separately. 



FIG. 12 : (Color online) (a) The y component my{x,0) and 
(b) the z component mz{x,0) of the local magnetization den¬ 
sity in the N for p = 7 r/ 2 , corresponding to the case where the 
magnetization in FI is perpendicular to that in F2. Other pa¬ 
rameters are the same as in Fig.fO] For comparison, ni}y(^^^(x) 
and m^y(^^^{x,0) are also plotted separately. 
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